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HIGHLIGHTS 


•  Online,  adaptive  RLS  and  TLS  estimates  of  battery  capacity  are  derived  in  closed  form. 

•  The  derivations  include  exact  weights  and  covariance  matrices  that  improve  accuracy. 

•  A  second  method  of  capacity  estimation  is  demonstrated  by  exploiting  battery  rests. 

•  An  optimal  fusion  algorithm  combines  independent  estimates  of  capacity  for  robustness. 
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In  this  paper,  the  second  of  a  series  on  battery  fuel  gauging,  we  present  an  approach  for  real  time 
capacity  estimation.  In  part  I  of  this  series,  we  presented  a  real  time  parameter  estimation  approach 
for  various  battery  equivalent  models.  The  proposed  capacity  estimation  scheme  has  the  following 
novel  features:  it  employes  total  least  squares  (TLS)  estimation  in  order  to  account  for  uncertainties  in 
both  model  and  the  observations  in  capacity  estimation.  The  TLS  method  can  adaptively  track  changes 
in  battery  capacity.  We  propose  a  second  approach  to  estimate  battery  capacity  by  exploiting 
rest  states  in  the  battery.  This  approach  is  devised  to  minimize  the  effect  of  hysteresis  in  capacity 
estimation.  Finally,  we  propose  a  novel  approach  for  optimally  fusing  capacity  estimates  obtained 
through  different  methods.  Then,  the  proposed  algorithm  was  validated  using  hardware-in-the-loop 
(HIL)  data  collected  from  commercially  available  Li-ion  batteries.  The  proposed  approach  performs 
within  1%  or  better  accuracy  in  terms  of  capacity  estimation  based  on  both  simulated  as  well  as  HIL 
evaluations. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

The  state  of  charge  (SOC)  of  the  battery,  defined  as. 

_  Remaining  charge  in  Ampere  hours  (Ah) 

~~  Capacity  of  the  battery  in  Ampere  hours  (Ah)  ^ 

provides  vital  information  about  the  status  of  the  battery. 
Numerous  approaches  can  be  found  in  the  literature  for  tracking 
the  SOC  of  the  battery.  The  knowledge  of  SOC  and  battery 
capacity  are  used  to  estimate  the  time  to  shut  down  (TTS)  or 
time  to  fully  charge  (TTF)  the  battery.  It  is  understood  that  the 
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battery  capacity  varies  with  temperature  and  that  it  fades  over 
time  depending  on  usage  patterns  and  age.  Accurate  tracking  of 
battery  capacity  is  a  critical  element  of  battery  fuel  gauging; 
however,  it  has  received  relatively  little  attention  in  the 
literature. 

Some  important  contributions  to  battery  capacity  estimation 
were  presented  in  Refs.  [1—3].  In  Ref.  [1]  a  method  of  estimating 
battery  capacity  based  on  SOC  lookup  of  a  resting  battery  and  the 
amount  of  transferred  charge  was  reported.  This  method  neglects 
the  hysteresis  effect  [4]  and  assumes  that  the  rested  battery  voltage 
represents  the  true  open  circuit  voltage  (OCV)  of  the  battery. 
Further,  the  proposed  approach  in  Ref.  [4]  is  not  designed  to  track 
capacity  drifts  over  time.  In  Ref.  [3  a  weighted  total  least  square 
(TLS)  estimation  scheme  was  proposed  for  online  battery  capacity 
estimation.  The  exact  value  of  the  weights  were  not  derived  in 
Ref.  3]  and  the  developed  TLS  method  was  not  directly  extendable 
to  adaptive  estimation.  In  Ref.  [2  ,  a  data  driven  approach  for 
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List  of  notations 

ftTLsk] 

TLS  capacity  estimation  error  variance,  (31) 

ws[k] 

process  noise,  (4) 

Q>att 

capacity  of  the  battery,  (2) 

ws[k] 

differential  SOC  error,  (12) 

Ch 

coulomb  counting  coefficient,  (6) 

—i< 

ws 

accumulative  of  process  noise  between  rest  at  k  and 

^LS 

LS  estimate  of  battery  capacity,  (22) 

the  subsequent  rest,  (40) 

CrLS 

RLS  estimate  of  battery  capacity,  (25) 

~,k 

^s,  OCV 

differential  SOC  error  based  on  OCV  lookup,  (44) 

Ctls 

TLS  estimate  of  the  battery  capacity,  (29) 

w; 

vector  of  differential  SOC  errors,  (16) 

d  c 

vector  of  SOC  differences,  (16) 

^S,  OCV 

vector  of  differential  OCV  based  SOC  errors,  (45) 

t\K 

ac,ocv 

vector  of  SOC  differences  between  two  rests,  (45) 

xc[k] 

fused  battery  capacity,  (51) 

d  t 

vector  of  Coulombs,  (16) 

xs[k\k} 

SOC  estimate,  (8) 

dsocv 

vector  of  Coulombs  between  two  rests,  (45) 

xs[k  +  1]  SOC,  (2) 

HK 

augmented  observation  matrix,  (26) 

^S,OCV  [k\ 

OCV  lookup  based  SOC  estimate,  (35) 

i\k] 

current  through  the  battery,  (3) 

^S.OCV  [k\ 

error  in  OCV  lookup  based  SOC  estimate,  (36) 

K, 

OCV  parameters:  I<0,  K\,  I<2,  I<3,  I<4,  I<5,  I<6,  I< 7,  (33) 

Zc[k] 

capacity  estimates  by  TLS,  RLS  or  OCV  method,  (52) 

rii[k\ 

noise  in  measured  current,  (3) 

Zj[k] 

measured  current,  (3) 

SOC  estimation  error  variance,  (15) 

ecbm  [k] 

capacity  estimation  error,  (58) 

QcM 

variance  of  capacity  drift,  (51) 

V  C 

charging  efficiency  (77  =  77J,  (7) 

RcM 

capacity  estimation  error  by  TLS,  RLS  or  OCV  method, 

V  d 

discharging  efficiency  (77  =  77^),  (7) 

(52) 

C X 

covariance  of  the  augmented  observation  matrix,  (27) 

ft  RLS  M 

RLS  capacity  estimation  error  variance,  (23) 

s* 

Ws 

covariance  of  differential  SOC  errors,  (20) 

concurrent  estimation  of  SOC,  parameters  an  battery  capacity  is 
presented. 

In  this  paper,  we  address  the  drawbacks  of  the  existing  ap¬ 
proaches  to  online  capacity  estimation  and  develop  an  approach 
with  the  following  novel  features: 

•  A  weighted  recursive  least  squares  (RLS)  estimate  of  capacity  with 
the  derivation  of  the  exact  weights.  Similarly  to  [3  ,  we  derive  a 
weighted  RLS  approach  for  online  capacity  estimation.  How¬ 
ever,  we  also  derive  exact  expressions  for  the  weights  based  on 
the  variances  and  covariances  across  times  of  the  updated  SOC 
tracking  and  the  current  measurement  error  standard 
deviation. 

•  A  novel  TLS  approach  for  real  time  tracking  of  battery  capacity.  We 
develop  a  novel  TLS  approach  that  gives  a  closed  form  expres¬ 
sion  for  capacity  estimation.  This  approach  can  be  used  for 
continuously  tracking  changes  in  the  battery  capacity. 


•  A  novel  method  for  adaptive  capacity  estimation  based  on  OCV 
lookup  of  a  rested  battery.  We  develop  a  TLS  approach  for  online 
tracking  of  battery  capacity  by  exploiting  battery  resting  in¬ 
stances  for  OCV  lookup  based  SOC  estimation. 

•  A  robust  approach  for  optimal  fusion  of  capacity  estimates  ob¬ 
tained  through  different  approaches.  We  develop  a  robust 
approach  for  optimally  and  adaptively  fusing  capacity  estimates 
obtained  through  two  different  approaches. 

In  Ref.  [5  ,  part  I  of  this  series,  we  presented  an  approach  for  real 
time  parameter  estimation  for  various  battery  equivalent  models, 
and  in  Ref.  [6  ,  part  III  of  this  series,  we  develop  an  SOC  tracking 
approach,  assuming  the  knowledge  of  parameter  and  capacity  es¬ 
timates.  Fig.  1  shows  a  block  diagram  of  the  BFG  and  the  steps 
followed  to  obtain  it.  We  briefly  describe  each  of  the  blocks  below. 

•  Offline  data  collection.  The  understanding  of  the  behaviors  of  a 
battery  is  achieved  through  a  combination  of  physical, 


Fig.  1.  Elements  of  a  BFG.  The  focus  of  the  present  paper,  capacity  estimation,  is  shaded  in  grey.  The  dashed  lines  indicate  that  the  CF/PF  models  are  not  used  in  tracking  SOH  and 
computing  RUL  in  the  present  paper  series. 
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electrochemical  model  based  explanations  and  through  data. 
For  example,  a  particular  battery  cell  having  a  certain  maximum 
voltage  and  the  fact  that  there  exists  a  monotonic  relationship 
between  the  OCV  and  SOC  of  the  battery  can  be  explained  based 
on  electrochemical  understanding;  the  modeling  of  OCV-SOC 
relationship  as  a  Combined+2  model  [7]  is  achieved  through 
data.  As  far  as  the  BFG  algorithm  development  is  concerned,  the 
data  collection  plays  an  important  part  in  understanding 
different  models,  such  as,  OCV-SOC,  CF  and  PF  models,  that  are 
required  for  accurate  fuel  gauging.  It  must  be  noted  that  data 
collection  and  battery  modeling  go  hand  in  hand  until  adequate 
models  are  defined  and  the  parameters  estimated. 

Battery  modeling.  There  are  different  paradigms  in  batter 
modeling  as  it  applies  to  fuel  gauging,  such  as  electrical  equiv¬ 
alent  circuit  models  (ECM)  based  approaches  5,6,8  ,  electro¬ 
chemical  modeling  based  approaches  [9]  and  data  driven 
approaches  [10-12  .  Our  BFG  approach  is  based  on  ECM,  hence, 
our  discussions  and  descriptions  are  related  to  ECM  modeling 
based  BFGs  where  the  following  five  aspects  of  the  battery 
behavior  are  studied  and  models  developed. 

1.  OCV.  The  characteristics  of  OCV-SOC  relationship  were 
extensively  studied  in  the  literature  (see  Ref.  [7]  and  the 
references  therein).  It  is  understood  that  the  OCV-SOC  rela¬ 
tionship  is  stable  over  temperature  changes  and  aging.  It 
must  be  noted  that  even  though  there  exists  a  clearly  un¬ 
derstood  OCV-SOC  characteristic  for  which  the  parameters 
can  be  accurately  estimated,  the  SOC  cannot  be  simply 
computed  based  on  the  measured  voltage;  the  dynamic 
behavior  of  the  battery  has  to  be  accounted  for  before 
computing  SOC  based  on  measured  voltage. 

2.  Battery  capacity.  Battery  capacity  denotes  the  total  amount  of 
Coulombs  that  can  be  discharged  from  the  battery.  Battery 
capacity  is  a  key  battery  parameter  in  computing  the  state  of 
charge  online,  through  Coulomb  counting. 

3.  Dynamic  behavior  of  the  battery.  The  measured  terminal 
voltage  of  the  battery  is  not  always  equal  to  the  OCV.  The 
voltage  drop  between  the  OCV  and  the  battery  terminal  is 
modeled  through  ECM  (also  see  Fig.  4). 

4.  Capacity  fade  (CF).  Battery  capacity  is  known  to  degrade  as  a 
result  of  aging,  temperature  and  DoD  (frequent  DoD  values 
closer  to  both  0  and  1  negatively  affect  capacity  13]).  Unlike 
OCV  modeling,  which  requires  data  from  just  one  char¬ 
ge-discharge  cycle,  CF  modeling  requires  extensive  data, 
over  several  hundreds  of  cycles,  for  analysis;  the  CF  modeling 
requires  continuous  data  from  new-to-dead  battery  under 
various  temperature,  aging  and  DoD  scenarios. 

5.  Power  fade  (PF).  Similar  to  CF,  the  ECM  parameters  change  as  a 
result  of  temperature,  aging  and  depth  of  discharge  (DoD) 
patterns  a  battery  undergoes.  Especially,  the  internal  re¬ 
sistances  increase  as  a  result  of  aging  and  this  limits  the 
amount  of  power  that  can  be  drawn  from  the  battery; 
resulting  in  power  fade.  In  addition,  it  must  be  noted  that  the 
available  power  varies  instantaneously  as  a  result  of  tem¬ 
perature  changes  and  this  adds  an  extra  layer  of  difficulty  in 
modeling. 

Offline  parameter  estimation. 

-  OCV  model  parameters.  For  OCV  parameter  estimation,  slow 
discharge  data  is  required  for  an  entire  discharge-charge 
cycle.  Several  models,  such  as  linear,  polynomial,  exponential 
and  combinations  of  these  have  been  proposed  in  the  liter¬ 
ature.  We  present  a  summary  of  the  existing  OCV  modeling 
approaches  in  Ref.  [7]. 

-  CF  parameters.  There  has  been  significant  interest  in  devel¬ 
oping  a  capacity  fade  model.  Several  data  analyses  and  at¬ 
tempts  at  creating  a  CF  model  can  be  found  in  Refs.  13-21  . 


Cycle  number 

Fig.  2.  Discharge  capacity  over  500  cycles.  The  battery  is  discharged  at  1.5  A  rate. 

Slower  discharge  rates  are  used  for  ten  consecutive  cycles  in  every  100  cycles. 

-  PF  parameters.  Some  examples  of  PF  modeling  and  parameter 
estimation  are  found  in  Refs.  [22-25]. 

•  Online  parameter  estimation  (system  ID).  Given  the  OCV  param¬ 
eters  of  the  battery,  the  ECM  parameters  and  the  battery  ca¬ 
pacity  can  be  estimated  online.  Indeed,  there  are  some  recent 
works  [2,26,27]  that  attempt  to  estimate  the  OCV  parameters 
online  as  well. 

-  OCV.  In  this  paper  we  assume  the  knowledge  of  offline  esti¬ 
mated  OCV  parameters  that  are  used  in  the  estimation  of 
capacity  and  ECM  parameters.  In  Refs.  [2,26,27  ,  some  ap¬ 
proaches  are  developed  to  estimate  the  OCV  parameters 
online. 

-  Capacity.  The  focus  of  the  present  paper  is  online  capacity 
estimation.  We  develop  two  different  approaches  for  online 
capacity  estimation  and  present  an  optimal  fusion  technique 
that  results  in  a  robust  strategy  for  battery  capacity 
estimation. 

-  ECM  parameters.  We  develop  a  novel  strategy  for  the  esti¬ 
mation  of  ECM  parameters  in  Ref.  [5];  the  developed 
approach  is  significantly  novel  compared  to  the  existing 
approaches. 

•  Tracking. 

-  SOC.  A  clear  understanding  of  the  OCV-SOC  models  led  to  the 
development  of  SOC  tracking  as  a  recursive  Bayesian  esti¬ 
mation  problem  that  can  be  solved  through  non-linear  filters 
such  as  extended  Kalman  filter.  The  Coulomb  counting  model 
represents  the  instantaneous  changes  in  SOC  and  serves  as 
the  process  (plant)  model;  the  OCV  +  ECM  model  establishes 
a  relationship  between  the  measured  voltage  and  SOC, 
serving  as  the  observation  model  for  SOC  tracking. 

-  SOH  (PF  and  CF).  Unlike  SOC,  the  SOH  changes  slowly,  hence, 
there  is  no  need  to  track  the  SOH  at  the  same  frequency  as 
SOC.  Several  models  developed  in  Refs.  13-21  present  CF 
and  PF  in  terms  of  the  number  of  cycles.  Unlike  OCV-SOC 
relationship,  the  CF/PF-cycle  relationship  is  much  more 
complicated;  the  cycle  cannot  be  uniquely  defined;  the  ef¬ 
fects  of  temperature  exposure,  depth  of  discharge  (DoD) 
significantly  affect  the  model;  the  patterns  of  temperature 
exposure  and  DoD  are  not  uniquely  defined.  Hence,  online 
estimated  capacity  and  ECM  parameter  values  are  well  suited 
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Fig.  3.  Generic  description  of  OCV  lookup  based  SOC  estimates  and  the  OCV  lookup  error.  The  solid  line  in  blue  shows  the  true  OCV-SOC  characteristics.  The  dashed  lines  in  red 
describe  the  effect  of  hysteresis  when  the  battery  comes  to  rest  after  charging  or  discharging:  The  dashed  line  in  red  above  the  true  OCV-SOC  curve  (solid  line)  corresponds  to 
hysteresis  after  charging  and  the  dashed  line  in  red  below  the  true  OCV-SOC  curve  (solid  line)  corresponds  to  hysteresis  after  discharging.  For  a  given  rested  voltage  zv[k],  the 
corresponding  SOC  estimate  xs,ocv[/<]  can  be  obtained  by  inverting  (33).  Due  to  hysteresis  in  the  battery,  this  SOC  estimate  xS  Ocv[/<]  will  differ  from  the  actual  SOC  xs[k],  resulting  in  an 
OCV  lookup  error  Xs  ocv  [/<]  •  The  OCV  lookup  error  will  always  be  negative  during  discharging  and  positive  during  charging.(For  interpretation  of  the  references  to  colour  in  this  figure 
legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


for  battery  fuel  gauging.  It  must  be  noted  that  the  online 
estimated  capacity  and  ECM  parameters  can  be  improved 
with  the  help  of  CF/PF  models,  which  is  beyond  the  scope  of 
this  paper. 

•  Forecasting/Trending. 

-  TTS.  In  practical  battery  applications,  the  batteries  are  fitted 
with  protective  circuitry  that  switch  off  the  device  once  a 
particular  shut  down  voltage  is  reached.  Predicting  TTS 


requires  the  knowledge  of  the  ECM  parameters;  usually  the 
TTS  prediction  is  made  for  a  particular  load,  e.g.,  TTS  for  voice, 
TTS  for  video  etc. 

-  RUL  Similar  to  TTS,  the  RUL  predicts  the  remaining  cycles 
until  a  certain  SOH  is  reached. 

It  must  be  noted  that,  each  estimation/tracking/forecasting 
block  (module)  in  Fig.  1  must  be  performed  at  different  time  scales. 


Fig.  4.  Equivalent  circuit  model  (ECM)  of  the  battery.  The  open  circuit  voltage  (OCV)  is  denoted  by  Refs.  V0(s[/<]),  h[k\  denotes  the  hysteresis  of  the  battery  at  time  /<,  and  iN[k]  denotes 
the  current  through  the  resistor  RN.  The  relaxation  effect  becomes  zero  (voltage)  when  the  battery  is  rested  whereas  the  hysteresis  effect  remains  unaffected  by  resting  the  battery. 
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In  Table  1,  we  describe  the  time  scales  at  which  each  of  these 
modules  operate. 

It  is  important  to  distinguish  the  online  capacity  estimation 
method  presented  in  this  paper  (and  also  in  Refs.  1]  and  [3  ) 
against  any  CF  models  that  can  also  be  used  to  estimate  battery 
capacity  at  a  given  cycle.  For  an  example  of  CF  models,  let  us  refer  to 
Fig.  2,  which  shows  the  discharge  capacity  of  a  Li-ion  battery  over 
500  cycles.  A  simple  CF  model  can  be  obtained  by  fitting  the 
discharge  capacity  vs.  cycle  data  as  an  exponential  decay.  Now, 
given  the  cycle  number,  the  above  CF  model  can  be  used  for 
determining  the  capacity  of  the  battery.  The  discharge  experiment 
whose  results  are  described  in  Fig.  2  was  conducted  over  a  60  day 
period  during  which  the  battery  experienced  the  same  temperature 
and  charge/discharge  patterns.  The  suitability  of  the  above  model 
in  individual  battery  usage  situations  cannot  be  guaranteed.  An 
online  capacity  estimation  method  (such  as  the  one  developed  in 
this  paper)  does  not  depend  on  CF  models,  hence,  it  is  applicable  to 
any  battery  usage  conditions. 

The  rest  of  the  paper  is  organized  as  follows.  Our  procedure  for 
capacity  estimation  is  summarized  in  Section  2.  The  performance  of 
the  proposed  approach  is  demonstrated  through  experimental  re¬ 
sults  in  Section  3,  and  the  paper  is  concluded  in  Section  4. 


2.  Battery  capacity  estimation  and  fusion 


The  instantaneous  state  of  charge  (SOC)  of  the  battery  can  be 
written  as  the  following  process  model,  which  is  also  known  as  the 
Coulomb  counting  equation,  in  terms  of  the  measured  current  as 
follows: 


/<  +  ! 


k 


+ 


Tj  A 


3600Ct>att 


k 


+ 


k 


(2) 


where  xs[k]  <e  [0,1]  denotes  the  SOC  of  the  battery,  Cbatt  is  the  ca¬ 
pacity  of  the  battery  in  Ampere  hours  (Ah)  and  Z[[k]  is  the  measured 
current 


Zi[k\  =  i[k]  +  rxi[k] 


(3) 


as  =  C^aj 


where  the  Coulomb  counting  coefficient  is 


Ch  3600Cbatt 


Here,  tj  is  a  constant  that  depends  on  whether  the  battery  is  being 
charged  or  discharged,  i.e., 


r i  = 


Vc  ![/<]>  0 

vd  <  o 


and  A  is  the  (constant)  sampling  interval. 


2.2.  BFG-based  online  battery  capacity  estimation  using  RLS 

In  Ref.  6],  we  discuss  ways  to  estimate  the  SOC  based  on  both 
voltage  and  current  measurements.  Two  consecutive  SOC  values 
xs[k]  and  xs[k  +  1]  are  written  in  terms  of  their  estimates  as. 

xs[k\  =  xs[k\k\  +xs[k\k]  (8) 


x5  [/<  +  1  ]  =  xs[/<  +  l|/<  +  1]  +  xs[k  +  l|k  +  1] 


where  the  estimation  errors  xs[k\k]  and  xs[k  +  1 1/<  +  1]  have  zero 
mean  and  variances  Ps[k\k]  and  Ps[k  +  1  |/c  +  1],  respectively.  The 
covariance  between  two  consecutive  estimation  errors  is  [[28], 
problem  5-10] 

E{xs[k\k]xs[k  +  1  \k  +  1]}  -  (1  -  G[k  +  1  }H[k  +  l])Ps[k|k]  (10) 


where  G[k  +  1]  is  the  scalar  Kalman  gain  and  H[k  +  1]  is  the  scalar 
linearized  observation  model  at  time  k  +  1.  The  details  about  G 
[k  +  1]  and  H[k  +  1]  can  be  found  in  Ref.  [6]. 

Now,  let  us  rewrite  (2)  in  the  following  form. 


/<+! 


-Xs 


k 


?7A 

3600CbattZ‘ 


k 


+  Ws 


k 


(ii) 


which  is  corrupted  by  zero  mean  white  noise  n/[k]  with  standard 
deviation  (s.d.)  g\. 

The  process  noise  in  (2)  relates  to  the  measurement  noise  in  (3) 
as. 

ws[k]  = -chMii[k]  (4) 

and  is  zero  mean  with  s.d. 

Table  1 

Time  scales  of  various  state  and  parameter  estimation  in  BFG. 

Offline  Tiny  scale  Small  scale  Medium  Large 

(milliseconds)  (seconds/  scale  scale 

minutes)  (hours)  (cycles/ 

days) 


Offline 

OCV(SOC), 

parameter 

CF(DoD, 

estimation 

temperature, 

other), 

PF(DoD, 

temperature, 

other) 

Online 

ECM  Capacity 

parameter 

estimation 

parameters 

Tracking 

SOC 

SOH 

(PF  and  CF) 

Trending 

TTS 

RUL 

(forecasting) 


Substituting  (8)  and  (9)  in  (11)  yields. 


A  Xs 

k  +  1,k 

=  XS 

k  +  1 

k  +  1 

-Xs 

k 

k 

.  4ftA 

k 

+  Ws 

k 

—  r 

Gbatt 

(12) 


where 


V  h 


V 


3600 


(13) 


and  the  differential  error  given  by 

ws[k]  =  xs[k\k]  - xs[k  +  1  |fc  +  1]  +  ws[k] 

is  zero  mean  with  variance  (see  proof  in  Appendix  A) 


(14) 


E{ws[k]2}^s[k]=G[k+\]2S 


k  + 1 


(15) 


where  S[k  +  1]  is  the  innovation  covariance  of  the  Kalman  filter 
(The  details  about  the  Kalman  filter  for  SOC  tracking  can  be  found  in 
Ref.  [6]). 

By  considering  a  batch  of  Lc  of  samples,  (12)  can  be  written  in 
vector  from  as  follows. 

ds  =  Cba'udc  +  w's  (16) 

where  k  is  the  batch  number,  Lc  is  the  length  of  the  batch, 
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do 


[  Ax$  [kLc  —  Lc  +  1 ,  kLq  —  Lc\  AXs  [kLc  —  Lq  +  2 ,  kLq 
—Lc  +  1 . . .  Axs  [/cLc,  kLc  —  1]] 


dc  —  [^h^/cLc-Lc+l2f  [^c  —  1  Vh^-KLc-Lc+l^i  [K^c  —  Lc  +  2 


■■■VjAklMxLc]] 


w*  =  [w[kLc  -  Lc  +  1]  w[kLc  -  Lc  +  2]...w[/cLc]] 


and  w*  is  a  white  Gaussian  noise  vector  with  covariance 


(19) 


(20) 


which  is  an  Lc  x  Lc  diagonal  matrix  with  its  nth  diagonal  element 
given  by 


■>/c 

'w 


s/  nn 


G[kLc  -  Lc  +  n]2S 


kLc  —  Lr  +  n 


(21) 


Now,  the  LS  estimate  of  the  inverse  battery  capacity  is  given  by. 


cLS  =  KY  2 


K 

We 


W)  (S*. 


IK 


(22) 


and  the  variance  of  the  LS  inverse  capacity  estimate  is 


R 


RLS 


4  =  (  (dc)T(sK 


IK 


We 


(23) 


When  a  new  batch  of  {d*+1 ,  d*+1 }  pair  arrives,  the  LS  estimates 
can  be  recursively  updated  by. 


R 


RLS 


K  1 


XR 


RLS 


K 


+  (dc+1)  Kt1 


IK+1 


(24) 


— 1 

(  A”1 

- 1 

Crls 

K  1 

=  ft  RLS 

K  - 1-  1 

(^ftRLS 

K 

Crls 

K 

+  d 


r')r(tT1(1‘+1 


— 1 


where  Rrls[k]  is  the  Lc  x  Lc  information  matrix  for  capacity  esti¬ 
mation  and  A  is  the  fading  memory  constant. 

It  must  be  noted  that  dKc  in  (16)  is  constructed  from  measured 
current  values  that  are  known  to  be  noisy  whereas  the  LS  and  RLS 
estimation  methods  discussed  earlier  assume  that  dKc  is  perfectly 
known.  For  a  more  realistic  solution,  one  should  consider  the  un¬ 
certainty  in  dKc.  Next,  we  propose  an  approach  based  on  total  least 
squares  (TLS)  optimization  that  addresses  errors  in  dKc. 

2.2.  BFG-based  online  battery  capacity  estimation  using  adaptive 
TLS 

In  this  section,  we  develop  an  online  capacity  estimation 
method  based  on  TLS  which  assumes  uncertainty  in  both  d*  and  dKc 
in  (16).  For  more  details  about  the  TLS  method,  the  reader  is 
referred  to  [29]  and  the  references  therein. 

Let  us  construct  the  following  augmented  observation  matrix. 


H 


K 


d  s  d) 


(26) 


The  information  matrix  associated  with  the  augmented  obser¬ 
vation  matrix  is. 


Sjj,  =  (HK)rHK 


(27) 


Let  us  write  the  eigendecomposition  of  in  the  following  form. 
SJ,  =  VkAkVkT  (28) 

where 

-  Ak  is  a  diagonal  2x2  matrix  of  nonnegative  eigenvalues  ar¬ 
ranged  from  the  largest  to  the  smallest,  i.e.,  A*(l,l)  denotes  the 
largest  eigen  value  and  A/c(2,2)  denotes  the  smallest  eigenvalue. 

-  Each  column  of  the  2x2  matrix  Vk  =  [v'J,  v£]  has  the  corre¬ 
sponding  eigenvectors,  i.e.,  the  first  column  v\  is  the  eigenvector 
corresponding  to  the  largest  eigenvalue  and  the  second  column 

is  the  eigenvector  corresponding  to  the  smallest  eigenvalue. 

The  TLS  estimate  of  the  inverse  battery  capacity  is  then  given  by 
the  ratio  of  the  components  of  v^,  namely. 


C 


TLS 


K 


v|(l) 

V5(  2) 

S&(1,2) 


(29) 


5ft(l,l)-A*(2,2) 


where  v^i)  is  the  ith  element  of  and  Sjfj (i,j)  is  the  (ij)th  element 
of  Sfj.  The  derivation  of  29]  is  shown  in  Appendix  B. 

For  smoother  estimates,  the  information  matrix  in  (27)  can  be 
updated  with  a  fading  memory,  as  follows 


skh  =  XSK~ 


1  (H*/H 

H  + 


K 


Lc-  1 


(30) 


Now,  based  on  [30],  the  TLS  estimation  error  covariance  is 
(approximately) 


(25)  i?TLs[K 


1 


M 


(31) 


(zK)TSf,zK 

where  h-  is  the  ith  row  of  H\  M  is  the  number  of  rows  in  and 

V2O) 


v2  (2) 


—  1 


(32) 


2.3.  OCV-based  online  battery  capacity  estimation 


The  open  circuit  voltage  (OCV)  of  the  battery  can  be  written  as  a 
nonlinear  function  of  SOC  as. 


V0(xs[k\)  =  K0  + 


1<1  ,  I<2  ,  I<3  ,  I<4 


+ 


+ 


+ 


xs[L\  xs  [k\ 2  xs  [k] 3  xs[k]4 
+  /C6ln(xs [/<])  + 1<7 ln(l  -  xs[k}) 


+  ^5xs  [k] 


(33) 


where  the  coefficients  I<o ,  K\,  K2,  K 3,  K 4,  K 5,  I< 6  and  K7  can  be  esti¬ 
mated  offline  by  collecting  voltage  and  current  measurements  by 
slowly  charging  and  then  discharging  the  battery  (see  Ref.  [7  ).  A 
typical  OCV-SOC  characterization  curve  is  shown  in  Fig.  3  in  blue. 

To  illustrate  this  further,  consider  the  equivalent  circuit  model 
(ECM)  of  the  battery,  shown  in  Fig.  4.  The  OCV,  which  is  a  function 


1  This  is  an  approximate  LS  estimate  because  both  (17)  and  (18)  are  noisy.  An 
algorithm  based  on  TLS  is  described  in  the  next  subsection. 


This  numbering  is  again  specific  to  this  document. 
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of  SOC,  is  denoted  by  V0(s[/<])  where  s[k]  is  the  SOC  at  time  k.  The 
relaxation  effect  of  the  battery  is  modeled  as  a  series  of  RC  circuits; 
when  the  battery  is  fully  at  rest  the  relaxation  effect  become  zero, 
i.e.,  when  there  is  no  current  activity  through  the  battery  for  suf¬ 
ficient  time,  the  voltage  across  the  RC  circuits  becomes  zero.  “The 
hysteresis  in  Li-ion  batteries  is  generated  due  to  the  thermodynamic 
entropic  effects,  mechanical  stress,  and  microscopic  distortions  within 
the  active  electrode  materials  during  Lithium  insertion! extraction  [4]”. 
The  hysteresis  effect  does  not  vanish  when  resting  the  battery. 
When  the  battery  is  fully  rested,  the  measured  voltage  can  be 
written  as 


zv[k]  =  Vo  (s [/<] )  +  h[k\  (34) 

where  the  hysteresis  h[k]  corrupts  the  measured  OCV. 

The  OCV-SOC  characteristic  curve  can  be  used  to  get  a  measure 
of  SOC  whenever  the  battery  is  sufficiently  rested.  The  SOC  of  the 
battery  for  a  given  at-rest  terminal  voltage  (which  is  also  the  open 
circuit  voltage)  zv[1<\,  written  as 


x 


S,  OCV 


k 


fc 


1 


OCV-SOC 


(35) 


can  be  computed  using  the  OCV-SOC  characterization  by 
computing  the  inverse  of  (33).  There  are  several  methods  for 
computing  the  inverse  of  a  nonlinear  function,  such  as  Newton's 
method  and  binary  search  [31  .  We  also  refer  to  this  as  OCV  lookup 
based  SOC  estimate. 

The  SOC  estimate  in  (35)  is  corrupted  by  the  hysteresis  voltage 
as  follows 


XS  [k\  —  XS,OCV  [fc]  +  XS, OCV  [fc]  (36) 

where  the  OCV  lookup  error  xS  Ocv[fc]  is  caused  by  the  hysteresis 
effect  in  OCV  (see  Fig.  3  also).  It  must  be  noted  from  Fig.  3  that 
when  the  battery  comes  to  rest  at  time  k  after  a  discharging 
process,  the  OCV  lookup  error  xS  Ocv[fc]  will  always  be  negative. 
Similarly,  when  the  battery  comes  to  rest  at  time  k  after  a 
charging  process,  the  OCV  lookup  error  xs,ocv[fc]  will  always  be 
positive.  However,  the  magnitude  of  the  error  will  vary  with  the 
amount  of  hysteresis,  which  is  a  function  of  the  magnitude  of  the 
current  before  rest,  SOC  and  time. 

Now,  let  us  rewrite  (2  in  the  following  form 

Xs[fc  +  1]  =  Xs[k]  +  ChAZj[fc]  +  ws[fc] 
xs[k  +  2]  =  xs[k  +  1]  +  c^A/<+-|Zj  [fc  +  1]  +  ws[fc  +  1] 


(37) 

(38) 


Xs  [fc  T"  N\  —  x$  [fc  -T  N  —  1  ]  +  CfjA]^]\f_-[Zi  [fc  +  N  —  1  ]  -T  w$  [fc  -T  N  —  1  ] 

(39) 

By  adding  (37)  through  (39)  on  both  sides,  one  arrives  at  the 
following: 


/jk  —  f  1  /jk  _i_  y\/k 
us, ocv  batt  c,ocv  '  vvs, ocv 


where 


dk  = 

uc,  ocv 


3600 


OCV  XS,OCV 


fc  +  N 


X 


s,ocv 


fc 


(41) 


(42) 

(43) 


^S,OCV  —  XS,OCV  [fc]  —  XS,OCV  [fc  +  N\  +W^  (44) 

v - v - ' 

J< 

X 

/vd,  ocv 

It  must  be  noted  that,  regardless  of  the  fact  that  the  sign  of  OCV 
lookup  error  xs,0cv[fc]  is  biased  towards  the  battery  mode  <e  (char¬ 
ging,  discharging},  the  “differential  error”  xd  ocv  (defined  in  (44))  can 
either  be  positive  or  negative.  By  considering  a  large  number  of 
differential  errors,  we  assume  xkd  ocv  is  approximately  white. 

Let  us  assume  that  the  /cth  batch  of  differential  SOCs  d*ocv  is 

collected  between  a  first  set  of  rest  points  k  =  k\,  fc 2,...,  hc  and  a 
second  set  of  rest  points  k  +  N  =  1<\  +  N\,  k 2  +  Afe,...,  hc  +  Nlc, 
respectively,  and  given  by 


J/C 

as,ocv 


f "> — 1  JK  . 

Lbattac,ocv  '  ws,ocv 


,K 


where 


J/C 

”s,ocv 


dkl 

uS,OCV’ 


cr2 

uS,OCV’ 


dkic 

us,ocv 


(45) 


(46) 


J/C 

ac,ocv 


d/Cl  dk: 

UC,  OCV’  uc, 


OCV  7 


dl<lc 

uc,  ocv 


iT 


(47) 
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s,ocv 


w: 


ki 

s,  ocv 


wkl 

VVS,  OCV  5 
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l< 


Lc 


1  T 


s,  ocv 


(48) 


Now,  it  can  be  seen  that  (45)  is  in  the  same  form  as  (16)  with 
dj  ocv*  ^c,ocv»  an4  ocv  replacing  d^,  dKc,  and  w£,  respectively.  Hence, 
RLS  and  TLS  based  capacity  estimates  can  be  derived  for  OCV  based 
observations,  by  following  Sections  2.1  and  2.2  respectively.  Let  us 

denote  the  RLS  and  TLS  estimates  of  OCV  based  capacity  as  Cro[k] 
and  CT0  [k],  respectively. 

Remark:  It  must  be  noted  that,  given  Nr  number  of  rest  states 
interspersed  with  a  discharge,  one  can  make  Nr(Nr  -  l)/2  differ¬ 
ential  observations.  For  example,  for  Nr  =  4,  let  us  assume  the 
battery  is  in  rest  state  at  time  instances  t\,  t2, 1 3  and  fable  2  shows 
the  possible  rest  pairs  for  differential  observations. 

In  the  practical  usage  of  mobile  devices,  there  will  be  many  time 
epochs  of  rest  state  within  each  cycle  of  usage.  It  must  be  noted  that 
the  batch  length  Lc  can  be  time  varying. 


2.4.  Robust  capacity  estimation  through  fusion 


Xs  [fc  +  N] 


x< 


[i<]  + 


3600Ct>att 


k+N- 1 

+  Ws\j] 
j=l< 


(40) 


where  wk  is  zero  mean  with  standard  deviation  VNas. 

Assuming  that  the  battery  is  at  rest  at  time  k  and  k  +  N,  we  can 
write  (40)  as. 


In  this  section,  we  propose  an  approach  for  fusing  the  TLS  es¬ 
timates  of  capacity,  developed  in  Section  2.2  and  Section  2.3.  The 
TLS  approach  is  found  to  outperform  the  RLS  approach.  Hence,  the 
derivations  in  this  section  are  developed  for  TLS  based  capacity 
estimate.  The  exact  derivations  can  be  applied  for  fusing  RLS  based 
capacity  estimates  as  well. 

The  online  capacity  estimate  Ctls[k]  is  corrupted  by  the  uncer¬ 
tainty  in  dKc  (see  (16))  which  is  caused  by  the  errors  in  the  measured 
current  and  the  uncertainty  in  the  which  is  caused  by  the  errors 


956 


B.  Balasingam  et  al.  /  Journal  of  Power  Sources  269  (2014)  949-961 


in  the  SOC  tracking  algorithm.  Similarly,  the  OCV  based  capacity 
estimation  Cto[k]  is  corrupted  by  the  uncertainty  in  d£  ocv  (see  ^45  ) 
which  is  caused  by  the  errors  in  the  measured  current  and  the 
uncertainty  in  d*  that  is  caused  by  the  OCV  lookup  differential  error. 
We  assume  that  the  errors  in  online  capacity  estimates  et[/c]  and  the 
OCV  lookup  based  capacity  estimates  et0M  are  uncorrelated.  Based 
on  these  assumptions,  the  capacity  the  fusion  becomes  the  fusion  of 
two  independent  tracks  [[32],  Section  9.2]. 

First,  it  must  be  noted  that  we  estimated  the^inverse  capacities, 
i.e.,  the  capacity  estimates  CTLS[/c]  in  29)  and  CT0  [k]  in  Sections  2.2, 
2.3  respectively  are  estimates  of  1/Cbatt-  Correspondingly,  the 
respective  estimation  error  covariances,  Rtls[k]  in  (31)  and  Rto[k] 
correspond  to  the  inverse  capacity  estimates  as  well.  Based  on 
Taylor  series  expansion,  the  expected  value  of  the  TLS  based  dy¬ 
namic  capacity  estimate  and  the  corresponding  estimation  error 
variance  are  approximated  as  (see  Appendix  C  for  details) 


Table  3 

Computational  complexity  of  the  capacity  estimation. 

Method 

Complexity 

BFG  based  online  capacity  estimation  -  RLS  (computation  of  (25)) 
BFG  based  online  capacity  estimation  -  TLS  (computation  of  (29)) 
OCV  based  online  capacity  estimation 

Fusion  (computation  of  (56)  and  (57)) 

Total  complexity 

*{Lc) 

*(Lc) 

^  (Tr.ocv) 

constant 

& 


K 


(55) 


Now,  whenever  a  new  measurement  zc[k]  e{CTLS [/<'],  CT0 [/<  ']}  is 
received,  where  either  k'  =  k  or  k"  =  k,  the  fused  capacity  estimate  is 
obtained  as. 
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K  —  1 


K  —  1 
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Pc[k  -  1 1/<  -  1]  +  Qc[k  -  1 


x  (Zc[k\  -  Xc[k 


Pc  [x  —  1  |ac  —  1]  +  0c[K  —  1]  T~  Rc  M 

l|k-l]) 


(56) 


(50) 


where  CTLS  is  the  estimate  of  Cbatt  based  on  dynamic  data  and 
RtlsM  is  the  estimation  error  variance. 

By  following  the  same  procedure,  one  can  obtain  the  OCV  based 
capacity  estimates  Cjo[/c]  and  the  corresponding  estimation  error 
covariance  ftroM. 

Now,  let  us  assume  the  battery  capacity  to  be  a  random  variable 
that  undergoes  the  following  slowly  changing  Wiener  process. 


Xc[x  +  1]  —  xc[k]  +  WcM  (61) 

where  wc[/c]  is  assumed  to  be  zero  mean  white  Gaussian  noise  with 
variance  QcM- 

The  capacity  estimates  Ctls[k]  and  CT0  [k]  fit  the  following 
observation  model. 


Zc[k]  =Xc[k]  +  I1c[k]  (52) 

where 


where  xc[k  -  1  |/c  -  1]  is  the  previous  update  of  capacity  estimate 
and  Pc[k  -  1  \k  -  1  ]  is  the  previous  estimation  error  variance  which 
is  updated  as 

r  1  _  Rc[k\(Pc[k  —  1  |ic  —  1]  +  Qc[k  -  1]) 
c[  J  Pc[k-\\k-\]  +  Qc[k  -  !}  +  Rc{k]  {  J 

The  above  fusion  approach  can  be  similarly  used  for  fusing  RLS 
based  capacity  estimates. 

2.5.  Computational  complexity 

In  this  section,  we  discuss  the  complexity  involved  in  estimating 

the  battery  capacity.  Once  the  pairs  of  vectors  d£  and  d£  are  given 
the  BFG  based  online  capacity  estimates  are  obtained  by  computing 
(25)  and  (29)  for  RLS  and  TLS,  respectively.  The  complexity  of 

computing  the  SOC  values  dj  is  referred  to  the  SOC  tracking  block 

[6]  and  d£  can  be  computed  through  Lc  additions.  The  complexity  of 
computing  (25)  and  (29)  is  linear  in  Ref.  Lc;  The  lengths  of  the 

vectors  d£  icv  and  d£  icv  depend  on  the  number  of  rest  points  a 
battery  gets  during  a  cycle  and  usually  Lc>0 cv  Lc.  Table  3  sum¬ 
marizes  the  required  computational  complexity  for  capacity 
estimation. 


2c 


K 


(53) 


3.  Simulation  results 


k!  and  k  are  the  time  indices  of  latest  estimates  from  the  corre¬ 
sponding  algorithms  (TLS  and  TO,  respectively),  and  uc[k]  is 
assumed  to  be  a  zero  mean  white  noise  with  variance. 


In  this  section,  we  provide  performance  analysis  of  the  proposed 
capacity  estimation  method.  First,  we  assess  its  performance  using 
simulated  battery  data.  Then,  we  use  real  usage  data  collected 
through  hardware  in  the  loop  (HIL)  experiments  for  performance 
assessment. 


Table  2 

Possible  rest  pairs  for  OCV  based  capacity  estimation. 


First  rest  (/<) 

Second  rest  (k  +  N) 

i 

£2 

t-i 

£3 

h 

£4 

t-2 

£3 

t-2 

£4 

G 

£4 

3.2.  Performance  analysis  of  BFG-based  online  capacity  estimation 
method  using  simulated  data 

In  this  section  we  evaluate  the  performance  of  the  BFG-based 
online  capacity  estimation  methods  using  simulated  data.  We 
simulate  typical  battery  usage  data  by  following  the  details  pro¬ 
vided  in  Ref.  [5]:  First,  typical  battery  usage  profile,  consisting  of 
dynamic  discharge  and  charging  of  the  battery,  was  simulated  in 
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Fig.  5.  Simulated  data:  Details.  From  top  to  bottom:  Simulated  load  profile  i[k], 
computed  current  x^[k]  through  R i,  true  SOC  xs[k]  and  noisy  battery  terminal  voltage 
zv[k\. 


Fig.  7.  Simulated  data:  Capacity  estimation  performance  when  the  algorithm  is 
initialized  with  the  wrong  capacity.  The  algorithm  is  given  an  initial  capacity  of  1.4  Ah 
whereas  the  true  battery  capacity  is  1.5  Ah. 


terms  of  true  load  current  i[k\.  Then,  measured  current  was  simu¬ 
lated  based  on  (3)  by  adding  random  noise  with  s.d.  05  =  10~3  A. 
Finally,  the  battery  terminal  voltage  was  computed  according  to 
battery  equivalent  Model  3  discussed  in  Ref.  5  assuming  a  voltage 
measurement  error  s.d.  of  erv  =  1CT3  V  and  that  the  true  initial  SOC 
of  the  battery  is  xs[0]  =  1.  It  is  assumed  that  Cbatt  =  1.5  Ah, 
ric  =  0.998,  and  rid  =  1  during  the  simulation.  We  use  the  same  OCV 
parameter  values  and  dynamic  equivalent  model  parameter  values 
as  in  Ref.  [5]  for  the  simulation. 

In  Fig.  5,  we  show  the  simulated  battery  usage  data.  The  simu¬ 
lated  load  current  Zi[k]  is  shown  in  the  top  plot  and  the  simulated 
voltage  zv[k\  is  shown  in  the  middle  plot.  The  true  SOC  is  computed 
based  on  Coulomb  counting  (2)  by  assuming  the  knowledge  of  the 
noiseless  current  profile  i[k]  and  true  capacity  and  is  shown  in  the 
plot  at  the  bottom. 

The  simulated  voltage  zv[k\  and  current  information  z\k\  are  fed 
to  the  battery  fuel  gauge  presented  in  Ref.  [6  which  uses  the 
method  proposed  in  this  paper  for  capacity  estimation  and  our 


proposed  approach  in  Ref.  [5  for  parameter  estimation  for  adaptive 
SOC  tracking.  The  SOC  tracking  algorithm  uses  Lb  =  200  (20  s  of 
data)  for  parameter  estimation  and  Lc  =  30  for  capacity  estimation. 

Next,  we  consider  the  performance  of  the  algorithm  in  terms  of 
battery  capacity  estimation  by  defining  the  capacity  estimation 
error  as 


(58) 


where  Cbatt(=l-5Ah  in  this  case)  is  the  true  value  of  the  battery 

capacity  and  Cbatt[/<]  is  the  latest  estimated  battery  capacity  xc[k\k] 
at  time  k  of  the  mth  Monte-Carlo  run. 

Fig.  6  shows  the  performance  of  BFG  in  terms  of  capacity  esti¬ 
mation.  The  input  to  the  capacity  estimate,  CTLS  [#c] ,  is  denoted  as 
“Estimated  Capacity”.  The  estimated  capacity  (both  BFG  and  OCV 
based)  is  sent  to  the  fusion  algorithm.  Since  there  are  no  rest 
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Fig.  6.  Simulated  data:  Performance  comparison  of  capacity  estimation.  Top  plot 
shows  the  “Estimated  Capacity”  CTLS[/c]  against  time.  The  bottom  figure  shows  the 
percentage  error  in  capacity  estimation, 


Fig.  8.  Simulated  data:  Comparison  of  RLS  and  TLS  methods  for  capacity  estimation. 
Both  methods  use  LS2  method  for  parameter  estimation. 
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Fig.  9.  HIL  data:  Details  Top  plot  shows  the  terminal  voltage  of  the  battery  and  the 
second  plot  shows  the  load  current.  The  battery  model  is  Samsung  AB463651  and  the 
test  temperature  is  25C. 

periods  in  the  data,  all  the  capacity  estimates  are  BFG  based,  hence 
the  outputs  from  the  fusion  algorithm  is  labeled  as  “filtered  ca¬ 
pacity”  instead  of  “fused  capacity”;  fused  capacities  are  shown  in 
Section  3.3.  It  must  be  noted  that  xc[k\k]  is  updated  whenever  there 
is  a  new  capacity  measurement.  The  output  of  fusing  the  capacity 
estimates  xc[k\k\,  which  is  used  by  the  BFG  for  SOC  tracking,  is 
shown  as  solid  line.  In  other  words,  the  solid  line  at  any  give  time 
corresponds  to  Cbatt  used  by  the  SOC  tracking  algorithm.  In  this 
example,  the  algorithm  was  given  the  knowledge  of  true  capacity 
for  initialization.  In  the  next  example,  we  demonstrate  the  perfor¬ 
mance  with  a  wrong  capacity  initialization.  The  error  in  battery 
capacity  estimation,  eCbatt  [k]  is  plotted  in  the  bottom  of  Fig.  6. 

In  Fig.  7,  we  show  the  ability  of  the  algorithm  to  “catch  up”  if  the 
true  initial  capacity  is  not  known.  This  is  particularly  useful  in 
practical  applications  where  the  same  types  of  batteries  will  vary  in 
their  capacities. 

In  Fig.  8  we  compare  the  performance  of  RLS  and  TLS  methods 
for  capacity  estimation.  This  figure  shows  that  the  TLS  method  of 
capacity  estimation  is  better  than  the  RLS  method.  However,  both 
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Fig.  11.  Demonstration  of  capacity  estimation  through  fusion  of  BFG  and  OCV  based 
estimates.  The  top  two  plots  show  the  measured  voltage  and  current.  The  third  plot 
shows  the  estimated  capacities;  the  BFG  based  online  capacity  estimation  is  performed 
whenever  there  is  dynamic  current  through  the  battery;  the  OCV  based  online  capacity 
estimation  is  performed  based  on  consequent  battery  resting  periods.  Whenever  an 
online  capacity  estimate  becomes  available  (either  from  BFG  or  OCV)  it  is  used  for 
Fusion. 


Fig.  10.  HIL  data:  Battery  capacity  estimation.  The  true  battery  capacity  was  measured 
by  a  slow  discharge. 


methods  provide  capacity  estimates  that  converge  to  within  1  % 
accuracy  in  just  one  cycle  of  data. 

The  analysis  on  simulated  data  showed  that  the  TLS  estimate 
outperforms  the  RLS  estimate  of  battery  capacity.  As  a  result,  we 
will  consider  only  the  TLS  estimate  in  the  sequel  for  studying  the 
performance  analysis  using  real  data. 

3.2.  Performance  of  BFG-based  online  capacity  estimation  method 
analysis  using  HIL  data 

In  this  section,  the  BFG  based  online  capacity  estimation  algo¬ 
rithm  is  validated  through  experimental  data  collected  from 
commercially  available  Li-ion  battery  cells.  The  BFG  evaluation  was 
also  performed  at  a  certain  temperature  using  hardware-in-the- 
loop  (HIL)  test:  a  fully  charged  battery  was  placed  inside  a  ther¬ 
mal  chamber  which  was  brought  to  the  test  temperature  and  then 
the  battery  was  subjected  to  different  charge-discharge  current 
profiles  and  the  battery  terminal  voltages  were  collected  using  a 
data  acquisition  system.  In  this  paper,  we  limit  our  analysis  to 
constant  temperature  evaluations  only. 
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The  SOC  tracking  algorithm  developed  in  Ref.  [5  was  applied  on 
the  collected  voltage  and  current  data  at  each  temperature.  The 
OCV  model  parameters  {/<j}  for  these  experiments  were  obtained 
offline  by  using  a  similar  type  of  battery.  Based  on  the  configuration 
of  the  data  acquisition  system,  the  s.d.  of  the  voltage  measurement 
error  was  assumed  as  av  =  10-3  V  and  the  s.d.  of  the  current 
measurement  noise  was  assumed  to  be  =  1CT3 4  A.  The  charging 
and  discharging  efficiencies  were  observed  to  be  ric  =  0.998  and 
rid  =  1,  respectively. 

Fig.  9  shows  the  voltage  and  current  data  collected  during  HIL 
evaluations  using  a  Samsung  battery  with  serial  number  AB463651 
at  25  C. 

The  performance  of  BFG  in  terms  of  battery  capacity  estimation 
is  shown  in  Fig.  10.  The  true  battery  capacity  is  measured  through  a 
slow  discharge  test  at  the  same  temperature.  The  capacity  esti¬ 
mation  feature  of  the  BFG  is  shown  to  be  able  to  converge  to  the 
true  capacity  of  the  battery  within  3-4  hours  of  operation. 

3.3.  Performance  of  capacity  estimation  through  fusion  analysis 
using  HIL  data 

In  this  section,  we  demonstrate  the  capacity  fusion  method.  The 
first  two  subplots  of  Fig.  11  show  the  data  collected  from  a  Samsung 
Galaxy  (GS4)  battery  with  serial  number  L1G6LLAGS.  Starting  with 
an  empty  battery,  one  full  cycle  of  data  was  collected  by  first 
charging  the  battery  until  it  becomes  full  and  then  applying  dy¬ 
namic  loads  separated  by  two  rest  periods.  Then  the  battery  is 
applied  with  a  pulsed  current  of  240  mA  average  magnitude  until 
the  terminal  voltage  reaches  3  V.  The  same  cycle  of  data  is  dupli¬ 
cated  (not  actually  collected  from  the  battery)  in  order  to  allow 
more  rest  points  for  the  demonstration  of  the  fusion  algorithm.  The 
outputs  of  the  BFG  based  online  capacity  estimation  method  (TLS 
approach)  and  the  OCV  based  online  capacity  estimation  method 
are  fused  by  the  optimal  fusion  algorithm  developed  in  Section  2.4. 

The  plot  in  the  bottom  of  Fig.  11  shows  the  performance  of  the 
Fusion  method.  The  BFG  based  online  capacity  estimates  are 
available  when  there  is  dynamic  load  in  the  battery.  The  OCV  based 
online  capacity  estimates  are  available  whenever  the  battery  be¬ 
comes  rested.  Whenever  there  is  a  new  estimate  (either  from  the 
BFG  method  or  OCV  method),  the  Fusion  routine  is  executed.  It 
must  be  noted  that  the  OCV  based  estimates  are  significantly 
different;  the  effect  of  hysteresis  could  be  attributed  to  this.  Since 
the  same  data  is  duplicated  for  several  cycles,  there  is  no  change  in 
the  OCV  based  capacity  estimates.  In  typical  smart  phone  applica¬ 
tions,  a  phone  could  be  considered  to  be  in  rest  while  it  “pings”  the 
tower  without  any  significant  applications  running  for  few  mi¬ 
nutes.  Hence,  there  could  be  significantly  more  number  of  rest 
points  within  one  cycle  for  OCV  based  capacity  estimation.  The 
objective  here  is  to  demonstrate  the  two  types  of  capacity  estimates 
and  show  the  progression  of  the  optimal  Fusion  algorithm. 

4.  Conclusions 

This  paper  presented  the  following  novel  features  for  battery 
capacity  estimation  towards  the  advancement  of  battery  fuel 
gauging: 

(i)  A  weighted  recursive  least  squares  (RLS)  estimate  of  capacity 
with  the  derivation  of  exact  weights.  We  computed  the  exact 
formulation  of  the  weights  based  on  the  SOC  tracking  error 
covariance  and  current  measurement  error  standard 
deviation. 

(ii)  A  novel  TLS  approach  for  real  time  tracking  of  battery  capacity. 
The  TLS  estimate  is  derived  in  closed  form  and  can  be  used 


for  adaptive  estimation  by  updating  the  covariance  matrix 
with  a  fading  memory. 

(iii)  A  novel  method  for  adaptive  capacity  estimation  based  on  OCV 
lookup  of  a  rested  battery.  We  consider  the  source  of  OCV  lookup 
error  (hysteresis)  in  the  derivations  and  propose  an  approach 
for  adaptive  estimation  of  capacity  through  OCV  lookup. 

(iv)  A  robust  approach  for  optimal  fusion  of  capacity  estimates 
obtained  through  different  approaches.  Based  on  the  capacity 
estimates  and  the  estimation  error  covariances,  the  proposed 
approach  uses  a  Kalman  filter  for  adaptive,  optimal  fusion. 

The  applicability  of  the  proposed  approach  can  be  discussed 

under  four  different  scenarios: 

1.  Applicability  in  the  presence  of  temperature  changes.  Both  of  the 
methods  presented  in  this  paper  estimate  the  total  capacity  [7] 
of  the  battery.  In  Ref.  [7  we  define  the  total  capacity  as  the 
maximum  amount  of  Coulombs  that  can  be  discharged  from  a 
battery  at  cell  temperature  T.  Typical  capacity  estimation  takes 
from  several  minutes  to  a  few  hours  during  which  the  battery 
temperature  is  prone  to  change.  The  total  capacity  model  (in 
terms  of  temperature)  has  to  be  included  in  the  derivations  for 
improved  accuracy  during  temperature  changes.  We  limit  the 
focus  of  this  paper  to  fixed  temperature  only. 

2.  Applicability  in  the  presence  battery  aging.  Both  methods  of  ca¬ 
pacity  estimation  developed  in  this  paper,  BFG  and  OCV  based, 
assume  the  normalized  OCV  models  of  [7  .  The  normalized  OCV 
models  are  shown  to  be  stable  in  the  presence  of  aging.  Hence, 
the  proposed  capacity  estimation  procedure  will  perform  well 
regardless  of  the  age  of  the  battery. 

3.  Applicability  for  various  battery  chemistries.  The  OCV  model  is  the 
heart  and  soul  of  the  online  battery  capacity  estimation 
approach  developed  in  this  paper.  When  there  are  batteries  of 
different  chemistries  to  be  gauged,  there  is  an  uncertainty  in  the 
OCV  model.  In  Ref.  [33  ,  we  develop  a  battery  chemistry  adaptive 
BFG  where  the  BFG  keeps  a  library  of  possible  OCV  model  pa¬ 
rameters  and  uses  the  probabilistic  data  association  (PDA)  al¬ 
gorithm  [32  in  order  to  associate  the  correct  OCV  model  to  the 
battery  being  gauged;  we  demonstrate  the  ability  of  the  PDA 
algorithm  to  quickly  find  the  correct  OCV  model.  In  summary, 
the  proposed  online  capacity  estimation  methods  in  this  paper 
will  be  effective  when  they  are  applied  along  with  the  battery 
chemistry  adaptive  BFG  of  [33]. 

4.  Applicability  of  the  algorithm  for  varying  dynamic  load  conditions. 
Dynamic  loading  conditions  affect  both  capacity  estimation 
methods  differently.  Below,  we  provide  brief  details. 

•  Effect  of  dynamic  load  variations  on  BFG  based  online  capacity 
estimation:  In  Ref.  [5],  we  discuss  four  possible  ECMs  as 
suitable  for  four  different  dynamic  loading  conditions  of  the 
battery.  If  a  BFG  uses  a  single  ECM,  the  effect  of  different  dy¬ 
namic  loading  conditions  might  manifest  as  SOC  estimation 
error  due  to  model  mismatch;  this  might  affect  the  estimated 
capacity  as  well.  However,  it  must  be  noted  that  the  SOC 


3  It  must  be  noted  that  the  ECM  parameters  also  vary  with  temperature.  How¬ 
ever,  our  voltage  drop  model  based  derivation  in  Ref.  [5]  requires  only  a  few 
consecutive  samples  for  ECM  parameter  estimation.  Typically,  we  use  200  samples 
in  Ref.  [5]  for  parameter  estimation;  assuming  a  sampling  interval  of  100  ms,  just 
20  s  of  data  is  sufficient  for  accurate  parameter  estimation.  The  temperature  change 
in  Ref.  20  s  can  be  neglected  without  compromising  in  accuracy. 

4  In  Ref.  [7],  we  report  preliminary  results  that  show  that  the  OCV  parameters 
obtained  through  normalized  OCV  modeling  is  unchanged  over  aging  battery.  The 
aging  in  Ref.  [7]  occurred  through  unused  shelf  life  as  well  as  repeated  full-to- 
empty  cycling  using  the  same  load  profile.  The  effect  due  to  other  types  of  aging 
is  currently  being  investigated. 
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estimation  is  performed  at  a  tiny  timescale  (in  milliseconds, 
see  Fable  1 )  whereas  the  capacity  estimation  is  performed  at  a 
medium  timescale  (hours);  which  means  that  the  effect  of 
dynamic  load  variations  will  not  be  “felt”  that  much  at  the  BFG 
based  online  capacity  estimation  as  long  as  the  model 
mismatch  is  temporary  (In  Ref.  5],  we  propose  to  use  model  3, 
which  suits  the  battery  operations  during  major  portions  of 
the  time). 

•  Effect  of  dynamic  load  variations  on  OCV  based  online  ca¬ 
pacity  estimation:  The  hysteresis  is  a  function  of  load  current 
(among  other  factors  such  as  temperature,  SOC  and  the  path); 
hence  different  dynamic  loading  conditions  will  result  in 
different  OCV  lookup  error  values.  The  OCV  based  online  ca¬ 
pacity  estimation  approach  was  derived  for  this  kind  of  situ¬ 
ations;  the  adaptive  estimation  strategy  is  designed  to 
mitigate  the  effect  of  hysteresis  on  capacity  estimation. 
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Appendix  A.  Derivation  of  capacity  estimation  error 
covariance 


In  this  section,  we  derive  the  covariance  of  the  differential  error 
in  (14).  For  convenience,  the  differential  error  (14)  is  rewritten 
below. 

ws[k]  =  xs[k\k }  - xs[k  +  1 1 k  +  1]  +  ws[k } 

The  objective  is  to  compute  the  variance. 


E<  Wc 


k 


ws[k]Tl=R 


ws[k\ 


(A.1) 


Let  us  write  the  process  equation  (2)  in  the  following  forms: 


xs[k  +  1]  =  x5[k]  +  chAZi[k]  +  ws[k\ 


(A.2) 


Appendix  B.  Closed  form  derivation  of  TLS  capacity  estimate 


Let  us  write  the  2x2  matrix  A  =  Sfj  as 
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The  eigenvalues  of  A  can  be  written  as 
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where  X\  is  the  largest  eigenvalue  and  X2  is  the  smallest  eigenvalue. 
The  eigenvector  corresponding  to  A2  is 
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Appendix  C.  Transformation  of  the  inverse  capacity  estimates 


In  this  appendix,  we  present  an  approach  to  get  the  capacity 
estimate  and  the  estimation  error  variance  based  on  the  inverse 
estimate  and  the  inverse  estimation  error  variance.  Our  derivations 
are  based  on  34,35]. 

Let  us  assign  simple  variables  for  inverse  capacity  estimate  and 
the  error  variance,  i.e., 
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xs[k  +  1  \k  +  1]  =  xs[k\k\  +  chAZi[k]  +  G[k  +  \}v[k  +  1]  (A.3) 


where  xs[k  +  1  \k  +  1]  is  the  Kalman  filter  estimate  of  xs[k  +  1],  v 
[k  +  1]  is  the  filter  innovation  and  G[k  +  1]  is  the  Kalman  gain. 
The  difference  between  (A.2)  and  (A.3)  is 

xs[k  +  1]  =  xs[k]  +  ws[k\  -  G[k  +  1  ]v[k  +  1]  (A.4) 

which  can  be  rearranged  to  the  following  form 

xs[k\  -  xs[k  +  1]  +  ws [k\  =  ws[k\  =  G[k  +  1  ]v[k  +  1]  (A.5) 

Hence, 


G[/<  +  l]2S 


/<+ 1 


(A.6) 


Rx^£{(x-x0)2} 

Defining 


(C.3) 


y=m  =  l  (C.4) 

A 

our  objective  is  to  find  approximations  for  E{y]  and  F{(y  -  E{y})2}. 
Appendix  CA.  Expected  value  of  y 

The  second  order  Taylor  series  approximation  is  given  by 

y=f{x)  =f{x  o)  +/'(x0)(x-x0)  +lf"(x0)(x-x0)2  (C.5) 


where  S[k  +  1]  is  the  innovation  covariance. 


The  second  order  approximation  of  E{y]  is  given  by 
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E{y}  =  £{/(x o)}  +/'(x0)  (£{x}  -  £{x0})  +^"(*0 )f{(x  -  x0)2} 


J_  Rx 
xo+x30 

(C.6) 


Appendix  C.2.  Variance  of  the  expected  value  ofy 

Let  us  expand  f[x)  as  a  first  order  Taylor  series  around  the  true 
value,  xq. 


y  =f(x)  =  /(X0)  +f'(Xo)(X  -  *o) 

The  variance  of  y  is  given  by 

£{(y-£{y})2}=£{(f'(x0)(x-x0)): 


(C.7) 
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x 


0 


Now,  the  expected  values  of  the  capacity  estimate  and  its  esti¬ 
mation  error  variance  are  given  by 
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